Astronomy & Astrophysics manuscript no. 9795 


© ESO 2008 


October 15, 2008 





An alternative parameterisation for binary-lens 
caustic-crossing events 



A. Cassan 



OO 

O 

o 

(N 
•*-> 
o 

O 

in 



43 

Oh 
6 



(N 
> 

(N 

in 

00 

o 

OO 

o 



X 



Astronomisches Rechen-Institut (ARI), Zentrum fur Astronomie (ZAH), Heidelberg University, MonchhofstraGe 12-14, 69120 
Heidelberg, Germany 



Received ; accepted 



ABSTRACT 



Aims. Microlensing events are being discovered and alerted by the two survey teams OGLE and MOA at an increasing rate. Around 
ten percent of these events involve binary lenses. Such events potentially contain much information on the physical properties of the 
observed binary systems, which can then be used for e.g. statistical studies on binary objects in the Galactic disk or bulge. However, 
such events are usually not straightforward to study, because the model equations are strongly non-linear and there are many local 
minima that can fool the search for the best solution if the parameter space is not inspected with great care. In this work an alternative 
parameterisation for the binary lens fitting problem is proposed, in which the parameters involved are defined to represent as closely 
as possible the caustic-crossing features observed in most binary lens light curves. Furthermore, we work out an extension of the 
method in order to make use of the straight line fold caustic approximation, when the latter applies for both the caustic entry and exit. 
Methods. We introduce an alternative parameterisation in order to confine the exploration of the parameter space to regions where 
the models only involve caustic crossing at the dates seen in the light curve. 

Results. We find that the proposed parameterisation provides more robustness to the light curve fitting process, in particular in 
avoiding a code to get stuck in false minima. 

Key words. Gravitational microlensing - Extrasolar planets - Stars: binary - Methods: numerical - Methods: mathematical 



1. Introduction 



Since the original idea of iPaczvnskil {1986) to use the gravita- 
tional lensing effect to track the presence of dark objects in the 
Milky Way, more than three thousand microlensing events have 
been discovered, mostly in the direction of the Galactic centre. 
A microlensing event occurs when the light from a background 
bright star (or source) is bent by the gravitational field of an in- 
tervening massive body (or lens) passing close to the observer 
line-of-sight. The phenomenon produces an apparent transient 
brightening of the source flux, resulting in a typical bel l-shaped 
light curve. Curr ent surveys opera ted by the OGLE dUdalskil 
I2003I) and MOA dBond et al.ll200ll) collaborations provide lists 
of ongoing microlensing events by monitoring more than 10 8 
stars towards the Galactic bulge on a daily basis, while (up to 
the 2007 bulge season) collaborations like PLANET, RoboNet 
or yt/FUN use ground-based networks of telescopes to perform a 
round-the-clock follow-up of selected events. 

Since the information on the lens characteristics is con- 
tained in the lens mass distribution rather than in its luminos- 
ity, microlensing naturally opens a new window in the explo- 
rat ion of invisible ob jects wi thin the Galaxy. Indeed, shortly af- 
ter lPaczvnskH d!986l) 's work. lMao & Paczvnskil d 199 lb proposed 
microlensing as a promising tool to hunt for extrasolar planets, 
a few years before the first pla nets were discovered around a 
pulsar dWolszczan & Fraill 19921) by pu lse timing, and around a 
solar-like star (Ma yor & Ouelozlll995l) by radial velocity mea- 
surements. 

Searching for planets is today the leading motivation of 
Galactic microlensing observations. Amongst the ~ 300 exo- 
planets discovered so far, most of them using the radial veloc- 
ity technique, microlensing contributes seven planets. Four of 



them are giant planets of a few Jupiter masses: M OA 2003-BLG- 
53/OGLE 20 05-BLG-235Lb dBond et all 120041) OGL E 2005- 
BLG-071Lb dUdalski et alJ 120051: iDong et al.l 120081) . plus a 
slightly less massive Saturn/J upiter planetary sy stem analogue 
(OGLE 2007-BLG- 1 09Lb+c, [GaudLetalJ Hool, orbiting low- 
mass dwarf stars. A similar number of super-Ear th planets have 
been found as well: OG LE 2005-BLG-390Lb dBeaulieu et al.l 
120061: iKubas et al.ll2008l) . a 5.5 Earth-mass planet at ~ 2.6 AUs 
which makes it the first ever rocky /icy planet discovered around 
a main sequence host beyo nd the snow line, OGLE 2005 -BLG- 
169Lb dGould et alj|2006h . 13 times more massive than Earth, 
and recently MOA 2007-BLG- 192Lb dBennett et alj|2008l) , one 
of the least massive planet discovered so far with an estimated 
mass of only ~ 3M e , probably orbiting a very late M or brown 
dwarf. 

While microlensing has so far contributed little to the large 
number of extrasolar planet discoveries, it probes a region in 
the mass/orbit diagram - namely cool low-mass planets (cold 
Neptunes and super-Earths) at few AUs from their stars - where 
no other competitive technique has the required sensitivity. 
Furthermore, microlensing is sensitive to extremely low mass 
bodies, so that planets of mass below that of the Earth could be 
reached by a dedicated space mission dBennett & Rhidl2002l) . 

Most of the microlensing events detected do not show ev- 
idence of the presence of planets around the lenses. For these 
so-called single lens events, modelling is relatively easy and it is 
possible to use a large number of events to draw statistical con- 
clusions on the abundanc e of planets. Such studies have been 
carried out in particular bv lGaudi et al.l d2002l) for a selection of 
PLANET events in the period 1995-199~9 ]Tsapras et alJ d2003l) 
with a sample based on OGLE events between 1998-2000, and 
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ISnodgrass et al.l d2004l) using data form the OGLE 2002 cam- 
paign. Nevertheless, the nature of many anomalous events re- 
mains unclear because their modelling requires much more work 
and computational resources. Up to now, such events have led 
to separate studies in order to derive constraints on the pres- 
ence of planets around the lens, as e.g. performed for the very 
high magnification eve nt (i.e very sensitive to planetary signals) 
OGLE 2007-BLG-050 dBatista et al.ll2008l) . 

The steadily increasing number of anomalous microlens- 
ing events is thus posing a great challenge to modellers, for 
the proper treatment of an individual case requires one to 
probe many models, hence demanding very efficient modelling 
schemes based on robust mathematical and numerical methods. 
A large fraction of binary lens events become unmodelled each 
seaso n. Since these represent close to 10% dJaroszvnski et al.l 
|2004|) of the overall sample (close to a thousand events a year), 
a special effort to provide improved and more automated mod- 
elling strategies is needed. 

Events that exhibit caustic crossings are usually very reward- 
ing in terms of extracting stronger constraints on the physics of 
the binary lens system, such as the masses and physical sep- 
arations of the individu al components (e.g. iJaroszvnski et al.l 
120051; iKubas et"aT1 120051) . The first microlensing event show- 
ing observati onal evidence of a binary lens caustic crossing is 
OGLE No. 7 (lUdalskiet all 19941) . and the first planet discovered 
by microlensing M OA 2003-BLG-53/OGLE 2005-BLG-235Lb 
dBond et alj|2004l) also involves a caustic crossing, emphasing 
the importance of such binary lens events. 

Although the light curve features associated with caustic 
crossings can be very pronounced and easily recognised, the un- 
derlying best model can be hard to find since the parameter space 
often is afflicted by many local minima, even if the light curve 
is well covered. This problem is aggravated by using parameter- 
isations in which the trajectory is modelled by parameters that 
are not well related to the light curve features. This is for ex- 
ample the case in the widely used binary lens standard param- 
eterisation using the angle between the binary lens axis and the 
minimum distance between the trajectory and the origin of the 
coordinate system. A main reason for this is that the equations 
of lensing are strongly non-linear and changing a given trajec- 
tory parameter can have very unexpected consequences for the 
shape of the light curve. This is explained by the presence of ex- 
tended gravitational caustics when the lens is composed of more 
than one body. When the source crosses a caustic, its flux expe- 
riences a very steep magnification gradient. Hence, even a small 
change in the trajectory can cause a very strong change in flux. 
A desired characteristic of an alternative parameterisation would 
be to directly relate the model parameters to visible features in 
the light curve, such as the caustic crossing dates, whatever the 
source trajectory. All probed configurations in the minimisation 
process would then display caustic crossings at their observed 
position on the light curve. We aim here to provide an alterna- 
tive parameterisation of the binary lens caustic crossing model 
that works in this way. 

The paper is organised as follows: in Sec. [2] we introduce 
the basic principles of the method and describe the new parame- 
terisation, in particular the choice of a curvilinear abscissa along 
the caustic lines. In Sec. [3] we detail the implementation of the 
method, and discuss the possible fitting strategies and the in- 
clusion of more elaborate binary lens models within this frame- 
work. We also tackle some important aspects arising from the 
introduced set of parameters. In Sec. [5] we work out in detail a 
possible way to implement the so-called straight line fold caus- 
tic approximation within the presented framework, which allows 



under some assumptions extremely efficient computations. We 
summarise in Sec. [6] our main results, and discuss the broader 
implication of fitting schemes and strategies to provide more au- 
tomated modelling software in the context of future generations 
of robotic telescopes. 

2. The method 

2.1. Caustic curves 

The fundamental equation of Galactic microlensing is the lens 
equation, which links the angular position of a point-source ob- 
ject to the angular positions of its multiple images in units of 6%, 
the angular Einstein radius (lEinsteinlll936l) : 



4GM (D S -D L 
D S D L 



(1) 



with £>s, Z\ the observer-source and -lens dista nces respect ively, 
and M the lens total mass. In the formalism o f lwittl(ll990h used 
here, the binary lens equation is written using the complex nota- 
tion: 
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where we have chosen the coordinate system so that the more 
massive component of the lens is located at the origin, with the 
secondary body located on the left at a separation^ d. Finally, 
q < 1 denotes the binary lens mass ratio. 

Multiple lenses in general harbour caustic lines, which are 
imaginary lines on the plane of the sky where the magnification 
of a perfect point-source diverges. Mathematically, they draw 
closed curves where the Jacobian of the transformation z — > £ 
vanishes. In the following and for brevity, we will refer to the 
term "caustic" not only for the lines themselves, but also for the 
closed curves they form. In the case of binary lenses, caustics 
exist in three different topologies , usually referre d as close, in- 
termediate and wide binaries (see iDom inik 1999, for an explicit 
discussion about these limits for extreme mass ratios), which are 
illustrated in Fig. Q] In the close binary case (on the left in the 
figure), there exits a central caustic near the more massive body, 
as well as two secondary caustics^ above and below the binary 
symmetry axis. The intermediate binary configuration (in the 
middle) displays a single caustic (also called a resonant caus- 
tic when d ~ 1), while the wide binary (on the right) involves a 
central caustic plus an isolated single secondary caustic, located 
between the two lens co mponents. 

It has been shown by lErdl & Schneider! d 19931) in the math- 
ematical framework of catastrophe theory that there exist close 
analytical formulae which determine the bifurcation values be- 
tween these regimes, which only depend on the lens separation 
and mass ratio (d, q): 



0_ + qf 
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d = 



(W /3 r 

1 + q 



(3) 



(4) 



1 The term "separation" has to be understood as the projected dis- 
tance between the two lens components on the plane of the sky, in units 
of8 E . 

2 In the case of planetary mass ratios (q <k 1), secondary caustics are 
usually referred as "planetary" caustics. 
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where d c (q) and d w (q) are respectively delimiting the close- 
intermediate and intermediate-wide binary regimes for a given 
value of the mass ratio q. The resulting frontiers are plotted in 
Fig.ffl 




0.5 1.0 1.5 2.0 
Separation (d) 



Fig. 1. The three topologies of binary lenses: the close binary 
lens (on the left), which involves a central caustic plus two off- 
axis small secondary caustics (plotted here for d = 0.8 and q — 
1CT 2 ), the intermediate binary (in the middle, assuming d — 1 
and q = 10~ 2 ) with a single caustic, and the wide binary (on 
the right, d — 1.6 and q = 10 2 ) whith a central caustic and an 
isolated secondary caustic. The figure also shows the bifurcation 
functions d c {q) and d w (q) between the three regimes, which are 
respectively displayed from left to right as the two vertical black 
lines. When the mass ratio tends to zero, the extension of the 
intermediate domain also tends to zero. 



2.2. The binary lens parameterisation 

In the simplest case, modelling a binary-lens microlensing event 
requires seven parameters, which are classically the lens separa- 
tion d and mass ratio q, the time to of the closest approach of the 
source to the origin (here the primary body), u$ the correspond- 
ing impact parameter in units of Be, ?e the Einstein crossing time, 
a the angle of the trajectory with respect to the x-axis, and the 
source radius p* in units of 6%. The corresponding parameterisa- 
tion is illustrated in Fig. [2] 

While this straightforward parameterisation (hereafter re- 
ferred to as "classical") is very well suited to describe single lens 
or non-caustic crossing binary lens events, it suffers from being 
rather disconnected from observable features in caustic-crossing 
light curves. The main reason is that the corresponding parame- 
ters contain no information on where the source crosses the caus- 
tics. As a consequence, widely exploring the classical parameter 
space in caustic-crossing events poses an important computa- 
tional challenge. Firstly, because a large proportion of the tested 
models do not produce caustic crossings located at the observed 
position on the light curve, and only a fraction of the calculated 
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Re(0 

Fig. 2. A binary lens generates caustics, as displayed here (red 
lines) for a system of two masses located on the x-axis with sep- 
aration d = 1 and mass ratio q = 10~ 2 (intermediate binary 
lens). The more massive (or primary) body is located at the ori- 
gin of the coordinate system, while the lighter component is on 
the left. In the simplest case, the source centre trajectory is clas- 
sically described (in thin letters) by a straight line that passes at 
time ?o at the minimum distance uq from the coordinate system, 
and makes an angle a with respect to the x-axis. In the proposed 
alternative parameterisation (in bold face), the source motion is 
expressed by means of the caustic entry and exit dates f entry and 

fexit, Which OCCUr at positions gentry (Sentry) and <Texit(Sexit) On the 

caustics. The corresponding model parameters s ently and s e xit re- 
sult from a proper definition of a curvilinear abscissa along the 
caustic lines (cf. Sec. 12.31 1. 



models contribute efficiently in the minimising numerical pro- 
cedure to find the solution. Secondly, caustic-crossing events 
are known to e xhibit very uneven x 2 hyper-surfaces (see such 
a^ 2 map in e.g. iKubas et al.ll2~005l) . and usually harbour a large 
number of local minima which should be found and c onsidered 
as co mpetitive models to explain the light curve (e.g. iDominikl 
119991) . Furthermore, in standard minimising routines, parameters 
vary within bounded intervals or around starting values, which 
may cause the routines to become stuck in a false minima if one 
(or more) of the parameters is not correctly optimised. Hence, 
parameterisations oriented towards observed features in the light 
curve help to circumvent strong correlations and parameter de- 
generacies. Additionally, a set of different starting points should 
be considered to ensure the completeness of the parameter space 
exploration. 

One of the requirements when building the new parameteri- 
sation is that the number of parameters should remain the same 
as for the classical one, but with parameters which reflect as 
much as possib le the caustic crossin g features that are seen in 
the light curve. lAlbrow et al.l (1 19991) proposed a parameterisa- 
tion adapted to model a light curve where one caustic crossing 
is clearly identified in the light curve, and for which a straight 
line caustic approximation can be applied. In the general con- 
text discussed in this section, we keep to the more general case 
where finite source effects are calculated with no approximation, 
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because the new parameterisation only affects the geometry of 
the problem and not the way to compute the magnification (see 
Sec. 14. II for more details). Nevertheless, we examine in Sec. [5] 
the inclusion of the straight line caustic approxi mation in our 
framew ork, thus extending the work initiated by lAlbrow et al.l 
(119991) . 

We therefore introduce in place of the four parameters uq, a, 
?e and to, the dates of the caustic crossing at the entry and exit 
dates, f ently and f ex ; t , as well as the two corresponding curvilinear 
abscissa s e ntry and s ex \ t along the caustic, so that £ entty = £(s ently ) 
and ^ ex i t = £(s e xit) are me P omts m the source plane where the 
trajectory intersects the caustic (cf. Fig. 0. Within a suitable 
range for these four parameters, a theoretical light curve will dis- 
play caustic crossing features around the observed entry and exit 
dates. The geometric conversion between the two set of param- 
eters is straightforward, since £ ently and £ e xit define the direction 
of the source motion. 



2.3. Curvilinear abscissa along the caustic curve 

In order to define a proper curvilinear abscissa, ones needs 
to build an appropriate parameterisation of the caustic curves. 
A possible option is to follow the method initiated by 
lErdl & Schneider! dl993]), for whom the caustics are described 
in a polar coordinate system (r, (ft) in the source plane. One then 
needs to solve for cos(0) for a given value of r, hence parameter- 
ising the caustic line by r. A working im plementation of such a 
parameterisation is described in detail bv lDomirTi k (2007) in the 
context of image contouring. 

In this paper, we work out an alternative method, based on 
a result obtained by IWittl d 1990t) : for a point-source f exactly 
located on a caustic, its corresponding point-like image z on the 
critical curve in the lens plane satisfies the equation: 



1 



1 + q 



1 

— + 



1 



(z + dy 



(5) 



where </> e [0, 2n] is an ad hoc parameter, independent of the 
choice of the coordinate system. For a given binary lens (d, q) 
configuration, we can compute for every <p in the range 0-2^ 
the caustic points by writing Eq. (0 as a complex polynomial of 
fourth degree: 



z 4 +2dz 3 + (d 2 - 



2 dp'* d 2 p'* 
1 + q 1 + q 



(6) 



For a given <p, this equation has four distinct complex solutions, 
which means the caustics are entirely built when <p runs from 
to 2n. Moreover, since the caustics never intersect in the case of 
a binary lens (contrary to a triple lens for example), there exists 
a unique value of (p for a given £ on the caustic (otherwise one 
would have (((pi) = (((pi) with tp\ + (p2, which would mean 
an intersection of the caustic lines). In this respect, Eq. (0 is 
optimal in the sense that each solution of the equation maps one 
and only one point on the caustic lines, and the caustics are fully 
drawn for < <p < 2n. 

We now examine how the caustics are drawn with (p for each 
of the three topologies. Since the caustics are symmetric with 
respect to the x-axis, we restrict ourselves to 3(0 > 0. By doing 
so, only two of the four solutions of Eq. (0 are taken into ac- 
count. The caustic region where 3(£) < is built by symmetry 
afterwards. The corresponding caustics are drawn in Fig. [3] to 
which the reader is referred in the following. 



- In the intermediate binary case, <p — leads to three points, 
A, B and C. When <p increases to 2n, the caustics are simul- 
taneously drawn by two distinct branches: A— >C and C— >B. 

- In the wide binary case, = correspond to the four points 
A, B, C and D. When (p increases to 27r, the central caustic 
is built with the single branch C— >B, while at the same time 
the secondary caustic is built by the branch A— >D. 

- In the close binary case, (p — correspond to A, D and C. 
When (p increases to In, the central caustic is built with the 
single branch A— >D, while at the same time the single branch 
C— >C draws the full secondary upper caustic (so that no hor- 
izontal reflection is needed). The lower secondary caustic is 
then the symmetric image of the upper one. 

A first naive choice would be to simply take (p as a curvilinear 
abscissa. But, as shown in the upper panel on the left of Fig. [3] 
a regular sampling on (p (i.e. 6<p = est) results in a very irregular 
linear sampling of the caustics. Most of the points accumulate 
near the cusps, which makes this parameter not well-suited. To 
avoid this effect, we define a curvilinear abscissa s which sat- 
isfies the property: 6s = est => \((s + 6s) - £(s)\ = est'. Since 
at first order ((s + 6s) £(s) + (d(/ds)6s, the condition to be 
fulfilled reads: \d(/ds\ = K, with K a constant, so that: 



= Kx 



ds 
~di> 



(7) 



The abscissa s can finally be obtained as a function of (p by inte- 
grating the derivative ds/dep, choosing s(<ft) monotonic: 



s((P) 



r ds i r 4 



dC 
d(p 



dep + 5(0) , 



(8) 



where the normalisation constant K depends on the choice of the 
ranges for <p and s. This is discussed in detail in Sec. for the 
three topologies. A remaining technical point here concerns the 
calculation of the derivative d£/d(p. By differentiating Eq. (0 
with respect to (p and substituting Eq. (0, we first obtain: 

# = * +e *^I ; (9) 
d(p d(p d(p 

we then get the expression of dz/d(p by writing 

*e-* = * x ± e -* = -ie* t (10) 
d(p d(p dz 

which leads, after differentiating Eq. (0 with respect to z, to 



dz i (z + d) 2 + q z 2 



2 (z + d) 3 +qz- 



(z + d)z, 



(11) 



and so d£/d(p can be calculated analytically given a point £ on a 
caustic (i.e. a point zona critical curve). 

In Fig. |3] with the exception of the upper panel on the left 
commented on previously, we have used a regular sampling of s 
to draw the caustic lines, following the above prescriptions. One 
can notice the regular spacing of the points along the caustics, 
even for the small number of sampling points chosen here (130 
for a whole caustic). When (p runs between and 27r, it draws 
a caustic which may contain cusp singularities. Because of the 
nature of suc h a singular ity, the tangential vector to the curve 
vanishes (e.g. lWitt)ll990l) . and thus the derivative of £ satisfies 
at the cusps locations: (d(/d(p) cu; , p = 0, or according to Eq. (0, 
(ds/d(p) cu; , p = 0. The cusps can be seen in Fig.0at the location 
of the inflexion points of s(<ft). 
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Fig. 3. The two left columns show the caustic patterns in the three possible topologies: intermediate (upper panels, d = 1.1 and 
q - 0.1), wide (middle panels, central and secondary caustics, d — 2 and q = 0.1) and close (central and upper secondary caustic, 
lower panels, d = 0.5 and q - 0.1). The primary lens is at the origin of the coordinate system, with the less massive body on the 
left. In the upper panel on the left, the sampled points on the caustic are computed assuming a regular spacing for <p in equation 
Eq. (O, for a total of 130 points (in practice this number is larger and is adjusted to match the desired precision in the sampling). 
For comparison and using the same number of points, the next panel on the right shows the sampling of the caustics using a regular 
spacing for the curvilinear abscissa s described in the paper, which results in a regular sampling of the caustics. The other plots 
display the wide and close caustics, also using a regular sampling for s. The three plots in the right most column show the function 
s((p), which converts <p to s according to Eq. (0 and the choices made in Sec. [3] 



3. Implementation 

As seen in the previous section, the caustics are built by gath- 
ering different branches obtained by moving tp from to 2n. In 
our fitting strategy, we choose to independently parameterise the 
central and secondary caustics, and make the choice that a com- 
plete given caustic is fully drawn for < s < 2. 



The abscissa s is obtained by first computing the caustics 
with a uniform sampling in <p e [0, 2n], and then computing nu- 
merically s(tp) using Eq. (©. The corresponding functions are 
plotted in the rightmost column of Fig. [3] for each of the three 
topologies. As stated before, the cusps singularities can be seen 
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as the inflexion points of s((f>), because (d£/ds) CU! . v = 0. In prac- 
tice: 

- In the intermediate binary case, since two branches are 
needed to build the upper part of the caustic, we choose 
< (f> < 4n, assuming that < (p < 2n builds the branch 
A— >C and 2n < <p < An the branch C— >B. Since we have 

< s < 1 for half a caustic, we set s((f> — 0) = and 
s(<p = 4n) = 1 in Eq. dS). The lower part of the caustic 

1 < s < 2 is built by a horizontal reflection. 

- In the wide binary case, we separately treat the central and 
the secondary caustic. Since in both cases a single branch is 
needed to build the upper part of the caustic corresponding 
to < s < 1, we choose s(<p = 0) = and s(<p = 2n) = 1. 
The lower part of the caustic 1 < s < 2 is once again built by 
an horizontal reflection. 

- In the close binary case, we also separately parameterise the 
central and secondary caustics. The central caustic parame- 
terisation is obtained in the same way as for the wide binary 
central caustic. On the other hand, the upper secondary caus- 
tic is fully drawn with the single branch C— >C, so that we 
have s(<p = 0) = and s(<p = 2n) = 2. The lower secondary 
caustic is just its horizontal reflection. 

Computing and parameterising the caustics is the most time 
consuming part of the method. Hence an efficient way to use 
this parameterisation is to perform the fitting procedure on a pre- 
defined {d, q) grid. This strategy is also strongly indicated when 
searching for all possible local minima within a reasonable Ax 2 
from the global minimum, which can possibly lead to competi- 
tive solutions. For example. iKubas et alj d2005l) found two local 
minima of almost the same goodness of fit for OGLE 2002- 
BLG-069L, one associated with a close binary and a central 
caustic crossing, and the other involving a wide binary and a 
secondary caustic crossing. 

In practice, for a given (d, q) grid point we determine to 
which topology it corresponds (cf. Eq. ©) and choose between 
the central and secondary caustics. A regular sampling of the 
caustics is then performed, using a regular spacing for s. The pa- 
rameters f entry and f ex it are then chosen so that during the fitting 
process they stay close (cf. below) to the values guessed from the 
light curve inspection. 

4. Fitting strategy 

4. 1. Extension to elaborated binary lens models 

A first important fact about the presented parameterisation is that 
it is purely geometrical, and so does not affect the way the mag- 
nification is computed. The magnification can be calculated via 
classical methods, such as magnification maps dRattenbury et al.l 
20021 IWambs ganssl l 1 997b. contouring methods dDominikil2007t 
Gould & G aucherel 1 19971) . improved ray shooting techniques 
dDong et al.l 120061) or special approximations of the magni- 
fication function. Thus, finite source effects, including limb- 
darkening, can be used within this framework. 

Secondly, since only the positions on the caustic and their 
respective dates are fitted, one can consider a s ource trajecto ry 
affected by effects like parallax (see in particular lGouldl2004l) or 
xallarap (motion of t he source due to a massive but usual ly un- 
seen companion, e.g. iJiang et al.l 120041: iGhosh et al. 20041). On e 
can also consider binary lens rotation (e.g. iDongetalJ 12008b . 
by computing and parameterising the (d, q) configurations at the 
entry and the exit dates, performing the appropriate coordinate 



Last, one can efficiently use the parameterisation to 
model events wi th poorly-sampled caustic-crossings (see e.g. 
iKains et aT1l2008l) : it prevents the code from converging towards 
a minimum that does not exhibit the caustic-crossing features 
at the right place, which is likely to happen with the classical 
parameterisation. Moreover it allows us to inspect different sce- 
narios with caustic-crossings occurring at differen t possible es- 
timat ed dates, such as for MOA 2007-BLG-197 dCassan et al.l 
120081) where the caustic entry date cannot be directly inferred 
from the collected data. 

4.2. Fitting using standards minimising algorithms 

When the parameterisation is used with minimising algorithms 
like the downhill simplex or conjugate gradient methods, which 
are well-suited to explore a small region around a local mini- 
mum, we choose the starting points for f entry and f ex it as estimated 
from the light curve inspection. The starting values for s entry and 
■Sexit (both evolving in the interval s e [0,2]), can be randomly 
chosen, or be set to values that allow us to explore trajectories 
of particular interest. A drawback of using these minimising al- 
gorithms is they may fail to locate all the minima, because they 
may converge to the nearest local minimum. This can be over- 
come by setting different sta rting points for s entrv and s eX rt- 

Genetic algorithms (e.g. ICharbonn eau 1995) are quite well- 
suited for use with the new parameterisation, because they com- 
bine the possibility of restricting the exploration of f ently and f e xit 
to intervals compatible with the positions of the caustic cross- 
ings on the light curve, along with a wide exploration of all kind 



of trajectories by performing a dense sampling of (s, 



entry ? 15 exit 

[0,2] x [0,2]. Because by essence a genetic algorithm does not 
converge to the exact minimum in the best found x 2 valley, it is 
recommended to further explore and refine the solution by using 
one of the method described above, or a Markov Chain Monte 
Carlo (MCMC, see below). 

The new parameterisation appears also very convenient for 
MCMC algorithms. The starting values for the parameters f entry 
and f ex i t can be estimated from the light curve as already men- 
tioned, while the one-cr intervals of variation around them can 
be adjusted to their estimated error. One can then perform a wide 
exploration of all possible combinations of caustic crossing posi- 
tions .sentry or s e xit, or refine the parameters around some already 
found local minima, thus choosing a small corresponding <x. We 
find in particular that the new parameterisation greatly increases 
the final mean acceptance rate. MCMC also has the great ad- 
vantage of providing the error bars on the best set of parameters 
without further calculations. 

Finally, an interesting property of the parameterisation is that 
the event time scale, fg, is computed rather than fitted as the usual 
case, so that 



fexit t. 



ts 



entry 



I gentry (fe3 



(12) 



rotation, and fitting e.g. for (d, 



entry* "exit? 'entry? 'exit? Gentry? ^exit 



for a straight-line and uniform source trajectory. This formula 
shows that there is a strong correlation between the distance be- 
tween the two caustic crossings and t^. Therefore it can happen 
that the minimising algorithm is not able to explore adequate 
values for fg, because for example the starting point was cho- 
sen too far from its final optimal value. It is then very likely to 
stop at a false minimum. Far from being a casual case, this is 
very likely to happen for binary lens caustic crossing events, for 
caustics show some similarity in their shape and magnification 
pattern over a wide range of binary separations and mass ratios 
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(see e.g. figures in lHarj2006tlChung et al.l2005tlDominikll999l) . 
which may draw a comparable magnification light curves, but for 
a range of values for t^. Using the new parameterisation, there is 
no such problem because fg is not constrained, which is one of 
its advantages. Finally, since the best models should be associ- 
ated with acceptable physical parameters, the values of fg found 
for the best model(s) can be used as a consistency check of their 
physical plausibility. 

5. The straight line caustic approximation 

5. 1. Assumptions and application 

The magnification in the vicinity of fold caustics has been 
stu died in detail by a number of authors (see i n particu- 
lar ISchneider & Weissl fl986t iGaudi & Pe tiers! l2002t iDominikl 
I2004T and one of the most important results in practice is that far 
enough from cusp singularities, the point-source magnification 
follows a generic behaviour. In particular, when a point-source 
crosses a caustic, a pair of critical images are created or merge at 
the critical line in the lens plane, depending on the source direc- 
tion. At first order, their magnification is the same, and is solely 
dependant on the point-source distance £ ± from the caustic, with 
A oc K I V?I, where the coefficient K is related to the local prop- 
erties of the lens equation. 

A first condition for a caustic to be treated as a straight line 
is that the local curvature radius of the causticQ should be much 
greater than the source radius. Moreover, the gradient of the 
magnification along the caustic should also be negligible, which 
is achieved when the source crosses the caustic far enough from 
a cusp. Within these assumptions, semi-analytical formulae can 
be derived which allow a very efficient computation of the time- 
dependant magnification during the caustic crossing, instead of 
the very demanding full finite-source integration. 

There are many possible applications of this approximation. 
In ord er to study microlensing effects in quasar lensing, for ex- 
ample, [%nehara| |200l|) makes use of a straight line fold caustic 
approximation to model a suspected caustic crossing affecting 
image "C" of quasar Q2237+03 05. In the context of stellar mi- 
crolensing, lAlbrow et al.l d 19991) have worked out an improved 
parameterisation for binary lenses involving one well-covered 
caustic crossing. In this framework, the new set of parameters 
is directly linked to the properties of the observed caustic cross- 
ing, such as the source radius caustic-crossing time, or the date 
at which it happens. This method was successfully applied to a 
binary lens event in t he Small Magellan ic Cloud microlensing, 
MACHO 98-SMC-l dAfonso et alfeoOOl) . Caustic crossings are 
also well-suited to study source stars' brightness profiles, be- 
cause different regio ns across thei r disks are differentially mag- 
nified: For example ICassan et al. I (12004 used the straight line 
caustic approximation to study the source limb-darkening for 
OGLE 2002-BLG-069 without requiring any knowledge of the 
full binary lens model. This is of great interest when the full 
light curve is not well covered, or when many different binary 
lens models fit the light curve. If the lens is a strong binary, with 
typical mass ratio q ~ 0.1 - 1, the caustics have a large spatial 
extension, which gives the method a high potential in practice. 

As mentioned in Sec. 12.21 we now examine in detail how 
the straight line caustic approximation can be introduced in 
our fram ework, hence provid ing an extension of the work ini- 
tiated by lAlbrow et al.l (1 19991) . In the following, we consider a 



given and fixed binary lens configuration (d, q). While outside 
of a caustic region a point-source is mapped onto three "non- 
critical" images, it is well known that a further pair of "critical" 
point-like images appear when the source crosses the c austic 
ilBourassa & Kantowskil 11975b ISchneider & Weissl Il986l) . The 
magnification for an individual image z issued from a point- 
source located at £ is given by 



A ps (z) 



1 



del I 



(13) 



where / = d(^,T])/d(x,y) is the Jacobian matrix associated with 
the lens equation Eq. (O, with z-x + iy and ( — % + ijj. 

Under the reasonable assumption that for many cases the 
source size is small with respect to the caustics, the three non- 
critical images are weakly affected by finite source effects dur- 
ing a caustic crossing. On the other hand, the magnification of 
the two critical images is strongly dependant on the finite exten- 
sion of the source. We thus divide the total magnification into 
two contributions, 



A(£p.)=A NC (f,p.)+A c (i,p.) 



(14) 



3 The local curvature of the caustic can be worked out analytically 
using in particular Eq. (9). 



where Anc stands for the three non-critical images total magni- 
fication and Ac for the two critical images total magnification, 
with p* the source radius in units of 6%. We detail in the follow- 
ing sub-sections how we can treat Anc and Ac in practice. 

5.2. The non-critical images magnification 

In previous studies making use of the straight line caustic ap- 
proximation, the total magnification of the non-critical images 
is approximated at first order by the Taylor expansion Anc - 
a + a) (t - t c ), where t c is the time of the caustic crossing (in the 
following, we use the subscript "c" for parameters computed on 
the caustic or critical lines, corresponding to the "entry" or "exit" 
used throughout the paper). Within this scheme, a and a> are fit- 
ted parameters. Th is choice is justified in alternative parameter- 
isations such as in lAlbrow et al.l (1 19991) . where only one caustic 
crossing is studied, or in ICassan et al.l (12004 where the caus- 
tic crossing is modelled alone, assuming no knowledge of the 
full binary lens model. In the parameterisation presented here, 
it would imply adding two more parameters per caustic cross- 
ing, for a total of four more parameters in addition to f en try> ?exit> 
Gentry, *exit> while the classical full binary lens parameterisation 
involves only one more parameter to treat finite source effects, 
p* (assuming a uniformly bright source). Moreover, it can hap- 
pen that for a caustic crossing located at the top of a microlens- 
ing light curve, the linear approximation of Anc is not precise 
enough because of a somehow stronger local curvature in the 
time-dependant magnification. 

Here we proceed in a different way: we choose to directly 
identify and track the three non-critical images during the whole 
caustic crossing, in order to compute their individual magnifi- 
cation. This can be done in the following way: at the position 
£ c where the source centre crosses the caustic line, we use the 
lens equation Eq. (O to locate the five images; the two critical 
images appear at the same position z c on the critical curve in the 
lens plane, with z c and ( c satisfying Eq. (0 as well. This provides 
a very robust numerical test to locate the two critical images - 
and therefore the three non-critical images at the caustic crossing 
date. A Newton-Raphson method can then be implemented to 
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"follow" the non-critical images and compute their point-source 
magnification 



A NC (<T,p« = 0) - J] A ps fc) 



(15) 



at any position during the caustic crossing. Alternatively and for 
more precision, one may replace the point-source approxima- 
tion by a Taylor expansion of the lens equation, and perform 
its integration over the source disk. Such an expansion of the 
m agnification up to Q(p 6 ) has been studied in the general case 
bv lPejcha & Hevrovskvl d2007l) . while a numerical implementa- 
tion up to the hexa decapol e for no n-critical images has been pre- 
sented in detail by iGouldl (120081) . The derived formulae do not 
require a further numerical integration, thus allowing a very ef- 
ficient computation of Anc(£>P* ^ 0), which is one of the main 
requirements about the met hod outlined in this section. In ad- 
dition to this, IGouldl d2008l) demonstrates that after few source 
radii from the caustic, the hexadecapole approximation is of very 
high precision. 

5.3. The critical images magnification 

For the critical images, the only relevant quantity in the straight 
line caustic crossing approximation is the perpendicular dis- 
tance from a point-source £ to the caustic. Following our choice 
of curvilinear abscissa along the caustic (cf. Sec. [3]), the inside re- 
gion delimited by the caustic curves is always on the right hand 
side of the tan gential vector, so that 4j_ > inside. As shown 
in particular bv lSchneider & Weissl d 1986b . the magnification of 
a point-like critical image in the vicinity of a caustic (and far 
enough from a cusp singularity) can be approximated to the first 
non-vanishing order by: 



A P C S (U > 0) * 



,1 



AH V2 (1-/1) 



1/2 



1 



A c s (^<0) = 0, 



(16) 
(17) 



where A = 1 - d^/dx and = dA/dy are evaluated at z = z c . 
Within these assumptions, the two critical point-like images 
are of equal magnification; the corresponding total finite source 
magnification is thus obtained by integrating 2 x A ps over the 
source disk. Assuming a uniform ly bright source disk, and fol- 
lowing ISchneider& Weissl (ll 9871) .' one obtains: 



A C (£,P») 



A 



WV2(1-A) 



1/2 



P, 



-1/2 j I l^-£'|sin</r c 



(18) 



where J is the sem i -analy tical function defined in Eq. (B.7) of 
ISchneider & Weissl d!987l) and f the centre of the source; if/ c is 
the angle between the source trajectory and the tangent to the 
caustic at critical point £ c . In practice, \j/ c = a — y c , where a is 
the trajectory angle of Fig.|2j and y c is the angle between the x- 
axis and the tangent to the caustic. Since the latter is collinear to 
(d£/d<p) :c , it can be evaluated analytically using Eq. d5). A linear 
limb-darkened source model can also be implemented very effi- 
ciently by adding one more p arameter per filter and computing 
further basis funct ions like J dAlbrow et al.ll 1 999t [Cassanll2005T 
ICassan et al]|2006l) . 

In practice, the function J can be easily tabulated to the re- 
quired precision, for it only depends on the ratio v - (±/p t . This 
is a key aspect of using the straight line caustic approximation, 
for one can replace a very demanding full finite-source compu- 
tation by a straightforward calculation of Eq. dT8~b , thus speeding 



up the computing time by many orders of magnitude. The tab- 
ulation of / can be restricted to a small interval, v = to few 
units (see below). 

5.4. The non-caustic crossing magnification 

Outside of the caustic crossing, the point-source or a finite- 
source approximation like the hexadecapole is used, while their 
validity for a given accuracy has to be estimated. In the sim- 
plest case of the point-source approximation, for example, the 
finite-source magnification of a critical image rapidly converges 
toward the classical point-source magnification: the absolute dif- 
ference between A ps and Ac is less than ~ 6 x 10~ 3 when the 
source centre is farther than three source radii from the caustic, 
and ~ 3 x 10~ 4 for ten radii. 

One way to ensure that the source crosses the caustic far 
enough from a cusp is to locate the curvilinear abscissa s c as 
done in Sec. 12.31 and to prevent the code from exploring values 
of s in the vicinity of s c . However if by chance the crossing oc- 
curs too close to a cusp, we expect a rather poor agreement of the 
magnification at the connection between the caustic crossing and 
non-caustic regions, which would naturally exclude the configu- 
ration obtained and thus minimise the effect of an inappropriate 
use of the approximation. 

Finally, within the straight line caustic approximation for 
both the caustic entry and exit, only one parameter remains to 
be fitted to account for finite source effects: the source radius p„, 
as in the classical parameterisation. 



6. Summary and prospects 

We have presented an alternative parameterisation of the caustic 
crossing binary lens problem, which has the main advantage of 
being directly linked to observed features in the light curve. We 
have explicitly discussed the underlying concepts and derived 
the corresponding equations, suggested a practical implementa- 
tion, as well as examined its implications in terms of numerical 
robustness. Furthermore, we have shown that the straight fold 
caustic approximation is particularly easy to include in the pre- 
sented framework, thus allowing a very high computation effi- 
ciency under some assumptions. 

Current survey and follow-up networks of telescopes gather 
microlensing data at an accelerating rate, making the modelling 
a rather challenging task. Moreover, observations now involve 
new generations of au tomated telescop es, such as the Rob oNet 
dBurgdorfetalJl2007l) and RoboNet-II dTsapras et alj|2008l) ini- 
tiatives, which will have in the near future the power to track 
many more microlensing events. In this respect, a more auto- 
mated modelling tool would be helpful to catch up with the flow 
of incoming data. Indeed, such a requirement for more auto- 
mated software also applies for telescopes based in isolated or 
remote locations, like Dome C in Antarctic a (see for example 
prospects underlined in lBeaulieu et alj|2005l) . where microlens- 
ing would benefit from a continuous monitoring of the Galactic 
bulge in addition to the site's climatic advantages. 

The automation of algorithms dedicated t o microlensing 
event prioritising has seen important efforts dDominik et al.l 
120071 ; ISnodgrass et al.l 120081) . but similar progresses in mod- 
elling tasks remain to be achieved on the mathematical and nu- 
merical sides to keep up with the increasing rate of detection of 
anomalous microlensing events, which is one of the motivations 
of this work. 
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